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Abstract 

We use the Bessel- inspired behavior of the structure function F2 at small x, obtained 
for a flat initial condition in the DGLAP evolution equations, with "frozen" and analytic 
modifications of the strong coupling constant to study precise combined Hl&ZEUS data 
for the structure function F2 published recently. 

1 Introduction 

A reasonable agreement between HERA data [I]-|S] and the next-to-leading-order (NLO) 
approximation of perturbative Quantum Chromodynamics (QCD) has been observed for 

> 2 GeV^ (see reviews in [4J and references therein) , which gives us a reason to believe 
that perturbative QCD is capable of describing the evolution of the structure function 
(SF) F2 and its derivatives down to very low values, where all the strong interactions 
are conventionally considered to be soft processes. 

A standard way to study the x behavior of quarks and gluons is to compare the data 
with the numerical solution to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equations [5] by fitting the parameters of x-profile of partons at some initial Qq and the 
QCD energy scale A [6l[7]. However, for the purpose of analyzing exclusively the small-x 
region, there is the alternative of doing a simpler analysis by using some of the existing 
analytical solutions to DGLAP equations in the small-x limit [3 [9]. 

The ZEUS and HI Collaborations have presented the new precise combined data [10] 
for the SF F2. The aim of this short paper is to compare the combined Hl&ZEUS data 
with the predictions obtained by using the so-called doubled asymptotic scaling (DAS) 
approach [9]. 

To improve the analysis at low values, it is important to consider the well-known 
infrared modifications of the strong coupling constant. We will use the "frozen" and 
analytic versions (see, |llj and references therein). 



2 Parton distributions and the structure function 

F2 

Here, for simplicity we consider only the leading order (LO) approximatior0. The struc- 
ture function F2 has the form 

F2(x,Q2) = efgix,Q'), fa{x,Q^) = f^{x,Q^) + f-{x,Q^), {a = q,g) (1) 
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where e = {J^i^D/f is an average charge squared. 

The small-x asymptotic expressions for parton densities look hke 
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where Ij^ (i^ = 0, 1) are the modified Bessel functions, 
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/3o' " ^ ' 27/3o' ^~ 27/3o 
denote singular and regular parts of the anomalous dimensions d+{n) and d-{n), respec- 
tively, in the limit n ^ lH. Here n is a variable in the Mellin space. 



3 "Frozen" and analytic coupling constants 

In order to improve an agreement at low values, the QCD coupling constant is mod- 
ified in the infrared region. We consider two modifications that effectively increase the 
argument of the coupling constant at low values (see [12] ) . 

In the first case, which is more phenomenological, we introduce freezing of the coupling 
constant by changing its argument — > + M^, where Mp is the p-meson mass (see 
[llj and discussions therein). Thus, in the formulae of Sec. 2 we have to carry out the 
following replacement: 

as{Q^) ^ afr(Q') = a,(Q2 + m^) (5) 

The second possibility follows the Shirkov-Solovtsov idea [T3] concerning the analytic- 
ity of the coupling constant that leads to additional power dependence of the latter. Then, 
in the formulae of the previous section the coupling constant as(Q^) should be replaced 
as follows: 

al^{Q')=as{Q')-^-J%- (6) 
Po Q - 

in the LO approximation and 

aan(Q^)=a.(Q^)-^^^ + ..., (7) 

in the NLO approximation. Here the the symbol . . . stands for the terms that provide 
negligible contributions when > 1 GeV |13j . 

Note here that the perturbative coupling constant as{Q'^) is different in the LO and 
NLO approximations. Indeed, from the renormalization group equation we can obtain 
the following equations for the coupling constant 

(8) 



in the LO approximation and 
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in the NLO approximation. Usually at the NLO level MS-scheme is used; therefore, below 
we apply A = Aj^. 



^ We denote the singular and regular parts of a given quantity kin) in the limit n — >■ 1 by k/{n — 1) and fc, 
respectively. 



Table 1: The results of LO and NLO fits to HI & ZEUS data fTO], with various lower cuts on 
Q^; in the fits the number of flavors / is fixed to 4. 
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Q2 > 5GeV^ 










LO 


623±0 055 


1.204±0.093 


0.437±0.022 


1.00 


LO&an. 


0.796±0.059 


1.103±0.095 


0.494±0.024 


0.85 


LO&fr. 


0.782±0.058 


1.110±0.094 


0.485±0.024 


0.82 


NLO 


-0.252±0.041 


1.335±0.100 


0.700±0.044 


1.05 


NLO&an. 


0.102±0.046 


1.029±0.106 


1.017±0.060 


0.74 


NLO&fr. 


-0.132±0.043 


1.219±0.102 


0.793±0.049 


0.86 


02 > 3 5GeV^ 










LO 


0.542±0.028 


1 089±0 055 


369±0 Oil 


1.73 


LO&an. 


0.758±0.031 


0.962±0.056 


0.433±0.013 


1.32 


LO&fr. 


0.775±0.031 


0.950±0.056 


0.432±0.013 


1.23 


NLO 


-0.310±0.021 


1.246±0.058 


556±0 023 


1.82 


NLO&an. 


0.116±0.024 


0.867±0.064 


0.909±0.330 


1.04 


NLO&fr. 


-0.135±0.022 


1.067±0.061 


0.678±0.026 


1.27 


> 2 5GcV^ 










LO 


0.526±0.023 


1.049±0.045 


0.352±0.009 


1.87 


LO&an. 


0.761±0.025 


0.919±0.046 


0.422±0.010 


1.38 


LO&fr. 


0.794±0.025 


0.900±0.047 


0.425±0.010 


1.30 


NLO 


-0.322±0.017 


1.212±0.048 


517±0 018 


2.00 


NLO&an. 


0.132±0.020 


0.825±0.053 


0.898±0.026 


1.09 


NLO&fr. 


-0.123±0.018 


1.016±0.051 


0.658±0.021 


1.31 


> 0.5GeV^ 










LO 


0.366±0.011 


1.052±0.016 


0.295±0.005 


5.74 


LO&an. 


0.665±0.012 


0.804±0.019 


0.356±0.006 


3.13 


LO&fr. 


0.874±0.012 


0.575±0.021 


0.368±0.006 


2.96 


NLO 


-0.443±0.008 


1.260±0.012 


0.387±0.010 


6.62 


NLO&an. 


0.121±0.008 


0.656±0.024 
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1.84 


NLO&fr. 


-0.071±0.007 


0.712±0.023 


0.529±0.011 


2.79 



4 Comparison with experimental data 

By using the results of the previous section we have analyzed Hl&ZEUS data for F2 
[10) . In order to keep the analysis as simple as possible, we fix / = 4 and as(M^) = 0.1168 
(i.e., A(^) = 284 MeV) in agreement with more recent ZEUS results given in til. 

As can be seen from Fig. 1 and Table 1, the twist- two approximation is reasonable 
for > 4 GeV^. At lower we observe that the fits in the cases with "frozen" and 
analytic strong coupling constants are very similar (see also dHHT]) and describe the data 
in the low region significantly better than the standard fit. Nevertheless, for < 
1.5 GeV^ there is stiU some disagreement with the data, which needs to be additionally 
studied. In particular, the Balitsky-Fadin-Kuraev-Lipatov (BFKL) resummation [15] 
may be important here [16]. It can be added in the generalized DAS approach according 
to the discussion in Ref. [17]. 

5 Conclusions 

We have studied the Q^-dependence of the structure function F2 at small-x values within 
the framework of perturbative QCD. Our twist- two results are well consistent with precise 
Hl&ZEUS data [10] in the region of > 4 GeV^, where perturbative theory is thought 
to be applicable. The usage of "frozen" and analytic modifications of the strong coupling 
constant, afr(Q^) and aan(Q^)) is seen to improve an agreement with experiment at low 



Q2 values, < 1.5 QeV^. 
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Figure 1: x dependence of F2{x,Q'^) in bins of Q^. The combined experimental data from HI 
and ZEUS Collaborations are compared with the NLO fits for > 0.5 GeV^ implemented 
with the standard (solid lines), frozen (dot-dashed lines), and analytic (dashed lines) versions 
of the strong coupling constant. 



